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Q\ , 1. Introduction 

O 

The similarities between self-consistent mean-field approaches to nuclear structure, such 
as Skyrme Hartree-Fock 0121, and Kohn-Sham density functional theory (DFT), which 
is widely applied to Coulomb many-body problems [3], are frequently noted. To go 
beyond the simple surface comparisons, however, there are many questions to address 
about applying DFT methods to nuclear structure and potential problems or challenges 
one will confront in trying to connect the two. For example, 

Ph ■ 

• How is Kohn-Sham DFT more than a mean-field calculation? Where are the 
^ . approximations? How do we include long-range effects in a Skyrme-like formalism? 

• What can you calculate in a DFT approach? Excited states? What about single- 
particle properties? 

• How does pairing work in DFT? Are higher-order contributions important? 

• How do we deal with broken symmetries (translation, rotation, . . . )? 

• Can we connect quantitatively to the free NN interaction? What about to chiral 
effective field theory (EFT) and/or low-momentum interactions? 

In this talk, we discuss how to make the nuclear connection to DFT systematically, 
by building a framework for addressing these questions. The idea is to use effective 
actions of composite operators to develop DFT [3] , with the Kohn-Sham orbitals arising 
through the inversion method used to effect the Legendre transformation to the effective 
action E] • The inversion method requires a hierarchy of approximations; here this is 
provided by an effective field theory expansion, which uses power counting to tell us what 
diagrams to include at each order |7. Ultimately, we seek nuclear energy functionals 
similar to conventional "phenomeno logical" approaches but model independent, with 
error estimates, and derivable from microscopic EFT approaches. 
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2. Effective Action Approach to EFT-Based Kohn-Sham DFT 

The dominant use of density functional theory (DFT) is in the description of interacting 
point electrons in the static potential of atomic nuclei. Applications include calculations 
of atoms, molecules, crystals, and surfaces [Sj. Discussions of DFT typically start with a 
proof by contradiction of the existence of the Hohenberg-Kohn (HK) energy functional 

E v [p] = F RK [p] + J d 3 a^(x)p(x) , (1) 

which is minimized at the exact ground-state energy by the exact ground-state density, 
and where Fhk[p] is universal. The Kohn-Sham procedure further posits a non- 
interacting system with the same density as the fully interacting system [3.. This leads 
to constructing E v using single-particle orbitals in a local potential i>ks( x )> 

N 

[-V 2 /2m + uks(x)]& = e t <pi => p(x) = £ |&(x)| 2 , (2) 

i=l 

which is determined self-consistently from E v [p\. The relative simplicity of solving 
Eq. (J2J for large finite systems is apparent while the entire approach is a win if there are 
reasonable approximations to Fhk[p] ; such as the local density approximation (LDA). 

We rely on a thermodynamic derivation of DFT, which uses the effective action 
formalism jH] to construct energy density functionals j3J |3] . The basic plan is to consider 
the partition function Z for the (finite) system of interest in the presence of external 
sources coupled to various quantities of interest (such as the density). We derive energy 
functionals of these quantities by Legendre transformations with respect to the sources. 
These sources probe, in a variational sense, configurations near the ground state. 

To derive conventional density functional theory, we consider an external source 
J{x) coupled to the density operator p(x) = tp^ (x)ip(x) in the partition function 

Z\J] = e~ w ^ ~ Tre-K S+J ® ~ / V[^]V[^\ e~^ c+J ^ , (3) 



for which we will construct a path integral representation with Lagrangian £ 8J. (Note: 
for convenience, we will take the inverse temperature /3 and the volume V equal to unity 
in the sequel.) The (time-dependent) density p{x) in the presence of J(x) is 

^ = J$ • M 

which we invert to find J[p] and then Legendre transform from J to p: 

5V[p\ 5T[p} 



T[p] = -W[J\ + Jjp with J(x) 



0. (5) 



Pgs(x) 



8p(x) Sp(x) 

For static p(x), T[p] is proportional to the DFT energy functional -Fhk' 

We still need a way to carry out the inversion; for this we rely on the inversion 
method of Fukuda et al. |1] . The idea is to expand the relevant quantities in a hierarchy, 

W[J, X) — Wq[J] + XWi[J] + X 2 W 2 [J] + ■■■ , (6) 

J[p, A] =J [p] + \Ji[p] + \ 2 J 2 {p] + • • • , (7) 

T[p, X] =T [p] + XTtlp) + A 2 r 2 [p] + ■ • • , (8) 
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treating p as order unity, and match order by order in A. Zeroth order is a noninteracting 
system with potential Jo(x): 

SW [J ] 



T [p] = -W [J 



d x Jo(x)p(x) 



p(x) 



(9) 



SJ (x) 

Because p appears only at zeroth order, it is given exactly from the non-interacting 
system according to Eq. 0. This is the Kohn-Sham system with the exact density! We 
diagonalize Wo [Jo] by introducing Kohn-Sham orbitals as in Eq. 0; it is the sum of 
Si's. Finally, we find Jo for the ground state by completing a self-consistency loop: 



j Q _> Wl _> Fl _> j x 



Wo 



J (x) 



5p(x) 



(10) 



i>0 i>0 

Note that even though solving for Kohn-Sham orbitals makes the approach look like 
mean field, the approximation to the energy and density is only in the truncation of 
Eq. (fTU|) at some order. Further, by using time-dependent sources we can generalize 
RPA to time-dependent DFT, to calculate properties of collective excitations. 

The hierarchy we have in mind is defined by EFT power counting. For example, 
the EFT for a dilute Fermi system with short-range interactions is defined by a general 
interaction as the sum of delta functions and their derivatives. In momentum space, 



(k|Kft|k'> = Co + ^C 2 (k 2 + k' 2 ) + C' 2 k ■ k' + • • • , 
which corresponds to £ e f t with general contact interactions (including 3+ bodies): 

2 vr T ' 16 



11} 



eft 



.d_ 

% dt + 2M 



1> 



^(^V) 2 + ^|[(^) t (^V 2 ^)+ h.c] 



+ ^(^v^) t -(^v^)- : ^(^W 

8 o 



+ 



(12) 



Dimensional analysis implies that ~ ^R 2l+l , Dn ~ ^i? 2l+4 , where R is the range 
of the interaction, which gives our hierarchy (e.g., if k F R 1). This £ ft will generate 
functionals like conventional Skyrme functionals (but with a different power counting). 
The conventional diagrammatic expansion of the propagator takes the form 



: 5 



+ 



+ 



+ 



: 3 



+ 



1 



+ 




£*fx,x';a;) 



with a non-local, state-dependent S*(x, x'; u). In contrast, in Kohn-Sham DFT we have 
a local Jq( x ) to all orders (i.e., not just at the Hartree level): 



= sr int [p] = ( f s P (x) \ 1 sr int [ P ] 

5p{x) 1 {5J (y)) SJo(y) 




+ 




+ 
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which introduces new Feynman rules for the "inverse density-density correlator" (double 
line) corresponding to [5p(x)/5J (y)] _1 - These inverse correlators appear in the energy, 



r«=00 + ^ + C<]>0+00=00+-" (13) 



but in this case the last two diagrams precisely cancel, in complete analogy to the 
cancellation of anomalous diagrams according to Kohn, Luttinger, and Ward |Hj. 

We can show how the full Green's function G is related to the Kohn-Sham Green's 
function G ks by adding a non-local source £(x, x') coupled to ip(x)ip^ (x') 0: 



z[J,Z] 



-W[J,£] 



(14) 



The Green's function follows as a functional derivative of W with respect to £ at constant 
J, which is in turn equal to a functional derivative of T with respect to £ at constant p. 
Since T can be decomposed as r from Eq. © plus everything else, 

r[p,£] = r [p,£] + r int [ P ,£] , (15) 

and Gks is the full Green's function for Fq, we can derive a Dyson-like equation: 



G(x,x') 



5W 



6T 



G ks (x,x') + G ks 



SG 



ks 



s P 



G 



ks 



(16) 



represented diagrammatically as: 



G 



G 



+ 



ks 



ks 



nG 



ks 




(17) 



We find a simple diagrammatic proof that the density from the exact Green's function 
(by taking a trace) equals the density of the non-interacting Kohn-Sham propagator: 




■O 



but other observables quantities of interest, such as single-particle properties, will 
generally differ and will need to be calculated using Eq. JT7J). 

To study further how well single-particle properties are reproduced, we include the 
kinetic energy density t(x) = (V^' ■ V?/>) in the DFT formulation for fermions in a 
trap. The meta-GGA DFT for the Coulomb problem uses r as an ingredient in the 
energy functionals, but applies a semi-classical expansion in terms of the density [3]. 
The result is a Kohn-Sham equation with a local potential only. In contrast, in Skyrme 
Hartree-Fock (SHF), p and r are treated as independent (N = Z nuclei here PQ): 



SHF 



[p, r, J] = J d 3 x{^r + h p 2 + ^3P 2+ ° + ^(3ti + 5t 2 )pr 



+ ^(9ti - 5t 2 )(Vp) 2 - ^W oP V ■ J + ^(h 



t 2 )J 2 



,(19) 
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(J is a spin-orbit density) and the orbital equations now have effective mass M*(x). 

To generate the analogous functional in DFT/EFT, add to the Lagrangian 
r](x) V^V^ and Legendre transform to an effective action of p and r: 



T\p,T]=W[J,r,]- Jj(x)p(x)- J 



7](x)t(x) . 



The inversion method results in two Kohn-Sham potentials, 

5T iQt [p,T] _ ^ cT int [p, r] 



Jo(x) 



and 



low 



5p(x) <5r(x) 
We diagonalize the quadratic part of the Lagrangian in Wo, 



d 4 x if) 



2M 



u(x) + Jo(x)- V-r/ (x)V 



and the Kohn-Sham equation becomes [with t>xs = u(x) — «/o( x )] 
r _ 1 



-V + uks(x)] if) a = e a ip 



(20) 



(21) 



(22) 



(23) 



L M*(xj 

with an effective mass l/2M*(x) = 1/2M — r] (x), just like in Skyrme HF 

To highlight the differences between an energy functional of p alone and one of p 
and r, we focus on the Hartree-Fock contributions from vertices with derivatives: 



-eft 



C2 [(#) t (^VV)+ h.c 1 ( ' 



16' 



+ 



The corresponding part of the energy density with spin degeneracy u 
LDA is 



2 (if}Vipy -(ipV<ip) + ... . (24) 
2 for Kohn-Sham 



£int[p] = • • • + 



Co 



2/3 



-8/3 



+ 



3C' 



6tt 5 



v 



2/3 



n 8/3 



+ ... (25) 



and with r is 11 



£mt[p,r] 



Vp) 2 ] + ^[pr-i(Vp) 2 ] + 



(26) 
-> t, 



4 8 4 

The pr functional with v — 4 takes the same form as the Skyrme functional with C ; 
(To get the spin-orbit part, one should couple a source to the spin-orbit density.) 

First, we note that power counting estimates work for all terms in the pr energy 
functional, as shown for two examples in Fig. ^ [the "pr" and "Vp" contributions are 
from Eq. 1)26)1 ]. The squares with error bars are estimated contributions to the energy 
per particle from terms in the energy functional at each order. The error bars indicate 
a "natural" range of coefficients between 1/2 and 2. The circles are the actual values. 
The hierarchy as well as the accuracy of the predictions is evident [TU] . 

We also note from these figures that the contributions that are new to the pr 
functional are very small. Thus, we are not surprised to find that predicted bulk 
observables are very similar (see Fig. |2J). The total binding energy and density are 
what Kohn-Sham is supposed to get right. What about the single-particle spectrum? 
The effective mass is unity in the p case and, as seen in Fig.^a), is reduced significantly 
in the pr case (the central values roughly cover the range in typical Skyrme interactions). 
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Figure 1. Contributions to the energy for two systems of trapped fermions |10j . 



The different effective mass is reffected in the single-particle spectra in Fig. Ofb), which 
can be understood qualitatively from the spectra calculated in a uniform system, 



7T 



£ k - £ k =-[iy- Ws + 2(1/ + i)o: 



k 2 - v 2 



P 



(27) 



I/-- • 2M 

Even though the density and energy are the same, the spectra are not. However, the pr 
result is closer to that of the full Green's function, as can be shown from Eq. (|17|) jlUj . 

An important ingredient of energy functionals for nuclei is pairing, which requires 
another generalization of Kohn-Sham DFT. To do so, we introduce an external source j 
coupled to the pair density, to break the U(l) phase symmetry associated with fermion 
number. The generating functional now becomes 



Z[JJ] 



-w[j,j] = I jju^^s e -S^ x [c + J{x)i>ti> ol +j{x){ii\i>l+^)] ^ ^8) 
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Figure 3. a) Effective mass and b) single-particle spectra for representative cases |1U| . 



and the fermion and pair densities p, <fi are found by functional derivatives with respect 
to J,j: 

6W[J,j] 



SJ(x) 



(29) 



(30) 



The effective action T[p, <p\ follows yet again by functional Legendre transformation: 

Y\p,<j>\ = W[J,j] - [ d A xJ{x)p{x) - [ d 4 xj{x)(f){x) , (31) 



and is proportional to the ground-state energy functional. The sources in turn are given 
by functionals derivatives with respect to p and 0, which are therefore determined in 
the ground state (where the sources are zero) by stationarity: 



5T[ P , 



Mx) 



5T[p, 



P=Pgs,' 







(32) 



P=Pgs,<f>=<t>g 



This is Hohenberg-Kohn DFT extended to pairing! 

Inversion is again carried out by order-by-order matching in an EFT expansion. 
Zeroth order is the Kohn-Sham system with potentials Jo( x ) an d Jo( x )) which yields 
the exact densities p(x) and 0(x). We introduce single-particle orbitals and solve 

h (x)-p j'o(x) \ / Ui(x) ) =E . ^( x ) 
jo(x) -h (x)+p J y Vi(-x) J 1 y Vi(-x) 

where Ui,Vi satisfy conventional orthonormality relations PQI2] and 



(33) 



/i (x) 



2M 



v x 



J (x) . 



(34) 



Thus it looks like conventional Hartree-Fock-Bogoliubov (HFB)! 

The diagrammatic expansion is the same as without pairing, but now uses Nambu- 
Gorkov matrix Green's functions. The Kohn-Sham self-consistency procedure is the 
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same as in Skyrme or RMF with pairing PQ I2]- I n terms of the orbitals, the fermion 
density is given by 



p(x) = 2^|^(x)| 2 , 

i 

and the (bare, unrenormalized) pair density is 

0bare(x) = [<(x)^(x) +Ui(x)v*(x)] = 2 ^ Uj (x)f j (x) . 

i i 

The chemical potential po is fixed by J p(x) = A and diagrams for T[p,< 
yields Kohn-Sham potentials 



Jn(x) 



5p(x) 



and j'o(x; 



P=Pgs 



<$r int [p, . 



(35) 

(36) 
-E\p,4] 

• (37) 



The pair density is divergent and requires renormalization. This problem is solved 
for a uniform system in Ref. ^1] by introducing a counterterm proportional to j 2 . In 
this uniform limit, is ultimately defined with a subtraction: 

" kc d 3 k 



(^) 3j0 \V(4-^) 2 +fi 4, 

which can be applied in a local density approximation (i.e., Thomas- Fermi), 



k c — >oo 



finite 



( x ) = 2^w i (x)v i (x) - j (x 



Mk c (x) 
2tt 2 



with 



2M 



+ -y KS (x) - p . 



(38) 



(39) 



However, convergence is very slow as the energy cutoff E c is increased, which is a 
practical problem for implementation in nuclei. Bulgac and Yu have shown how one 
can make different subtractions and greatly improve the convergence jT^j. For example, 
defining <p instead through 

" fcc d 3 k 



V 



■Jo 



k c — *-oo 



v^-^ + jo 2 4 - A*o 



finite 



(40) 



leads to much faster convergence, as seen here for the uniform system: 



10 
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--- Eq. (40) 
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Bulgac has found even faster convergence with other subtractions. 
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A question to be investigated is the role of higher-order contributions. In the weak- 
coupling limit, the pairing gap in a dilute Fermi gas with two spin states is reduced 
approximately fifty percent by induced interactions. The DFT/EFT framework is 
perfectly suited to explore the ramifications for finite nuclei, with a consistent treatment 
of interactions in the particle-particle and particle-hole channels. 



3. Problems and Challenges 

There are many problems and challenges facing the DFT/EFT program for nuclei. Here 
we give a personal laundry list of what we believe are priority issues. 

• Including long-range effects. Two classes of long-range phenomena are 
particularly relevant. On the left, we have long-range forces, such as pion exchange. 




On the right are diagrams leading to non-localities from near-on-shell particle-hole 
excitations. We will need appropriate expansions or resummations for both. 

• Developing gradient expansions. We need to consider semiclassical expansions 
used in Coulomb DFT and to revisit the density matrix expansion. Furthermore, 
there are gradient expansion techniques for one-loop effective actions, which might 
be adapted. The first targets are the beachball diagram for short-range forces and 
the Fock diagram for pion exchange. 

• Restoring broken symmetries. In Coulomb applications, a fixed external 
potential (nuclear Coulomb field) means that translational and rotational invariance 
are not issues, so they have not been studied via DFT. This is also true of particle- 
number violation in HFB. In the effective action framework, these correspond 
to problems from "zero modes," which can be dealt with using a Fadeev-Popov 
approach. Work in progress seeks an energy functional for the intrinsic density, 
using a one-dimensional toy model as a theoretical laboratory. 

• Developing auxiliary field Kohn-Sham theory. Auxiliary fields are like non- 
dynamical meson fields, which can be incorporated in the formalism using an 
appropriate saddlepoint evaluation. It may be profitable to revisit the large N 
expansion in this context. For pairing, the challenge is to separate particle- hole 
and particle-particle channels for short-range interactions. 

• Studying covariant DFT. Relativistic mean-field models also need to be cast into 
DFT form. We are exploring a controlled laboratory for finite systems using short- 
range interactions only. Special renormalizations allow for great simplifications. 
The connection to three-body forces in the nonrelativistic limit is being pursued, 
as well as covariant pairing and time-dependent Kohn-Sham theory. 
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Work in all of these areas is in progress. 

Perhaps the ultimate challenges for the DFT program for nuclear structure is to 
connect quantitatively to microscopic approaches. Past attempts to derive Skyrme 
parameters from NN interactions have resulted in qualitative or semi-quantitative 
agreement, but have always fallen short of useful quantitative predictions. In contrast, 
Kohn-Sham density functional theory for Coulomb systems is "ab initio" in the sense 
that the LDA part of the energy functional comes directly from numerical calculations 
(with a form based on analytic limits) of a uniform interacting electron gas. 

The successes of Coulomb DFT are often attributed to the fact that Hartree-Fock 
gives the dominate contribution, so that correlations are small corrections (even though 
the precision required — chemical accuracy - - is high!). In turn, it has been said 
that the same level of success is unlikely for nuclear systems since for conventional NN 
interactions, the correlations are larger than Hartree-Fock. However, these conventional 
potentials are not the only choice for many-body interactions. 

If we start with a chiral effective potential and then run a momentum cutoff down, 
we end up with a much softer potential. In the medium, as the cutoff is run toward 
the Fermi surface, the phase space for particle-particle scattering intermediate states in 
the medium becomes severely constrained. The apparent result is a more perturbative 
many-body system (at least in the particle-particle channel), so Hartree-Fock (with 
both two- and three-body low-momentum interactions) is a useful starting point 
Deriving constants for Skyrme-like energy functionals from this approach is an important 
problem for the near future. 
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